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ABSTRACT 

In low-coUisionality plasmas, anisotropic heat conduction due to a magnetic field leads 
to buoyancy instabilities for any nonzero temperature gradient. We study analogous in- 
stabilities in degenerate collisional plasmas, i.e., when the electron collision frequency 
is large compared to the electron cyclotron frequency. Although heat conduction is 
nearly isotropic in this limit, the small residual anisotropy ensures that collisional 
degenerate plasmas are also convectively unstable independent of the sign of the tem- 
perature gradient. We show that the range of wavelengths that are unstable is indepen- 
dent of the magnetic field strength, while the growth time increases with decreasing 
magnetic field strength. We discuss the application of these collisional buoyancy in- 
stabilities to white dwarfs and neutron stars. Magnetic tension and the low specific 
heat of a degenerate plasma significantly limit their effectiveness; the most promising 
venues for growth are in the liquid oceans of young, weakly magnetized neutron stars 
{B < 10^ G) and in the cores of young, high magnetic field white dwarfs {B 10^ G). 
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1 INTRODUCTION 

The presence of anisotropic thermal conduction in mag- 
netized low-collisio nality plasmas fundamentally changes 
the .Schwarzschildj (ligssl) criterion for convective stability 
l|Balbusll2000l . l200ll : lQuataertll2008l . hereafter Q08). Instead 
of convection being driven solely by an inwardly increasing 
entropy gradient, low-collisionality plasmas are convectively 
unstable for both outwardly and inwardly increasing tem- 
perature gradients (jBalbus .2000,, Q08). These instabilities 
have been dubbed the magnetothermal instability (MTI; 
dT/dr < 0) and the heat-flux driven buoyancy instability 
(HBI; dT/dr > 0). 

The key driver for both of these instabilities is the 
highly anisotropic heat flux established by even a dynam- 
ically weak magnetic field. In particular, if the electron cy- 
clotron frequency is large compared to the electron collision 
frequency, heat fiows almost entirely along magnetic field 
lines. However, in the presence of a magnetic field, the con- 
ductivity tensor always possesses some degree of anisotropy, 
even when collisions are rapid. Thus, any plasma in which 
electron conduction is energetically important may be sus- 
ceptible to the MTI and HBI. This motivates us to gen- 
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eralize previous analyses of conduction-mediated buoyancy 
instabilities to collisional plasmas. We are particularly inter- 
ested in the possible application of this physics to the highly 
conducting plasmas in neutron stars (NSs) and white dwarfs 
(WDs) and so we will allow for the possibility of degenerate 
plasmas. 

A reasonable measure of the coUisionality of a plasma 
is given by the ratio of the electron cyclotron frequency to 
the electron collision frequency, UgT, where cug = eB/^nieC 
is the Larmor frequency of electrons with Lorentz factor 7, 
B is the magnetic field strength, and r is the time between 
collisions. To motivate the parameter regime of interest in 
this paper, we estimate r as the inverse of the electron- 
ion collision rate, Uei = UiaeiVe, where rii is the number 
density of ions, a^i is the Coulomb cross section, and Ve is 
the (thermal or Fermi) velocity of the electrons. For the de- 
generate plasmas of interest in WDs and NSs, the kinetic 
energy equals the Fermi energy Ep. Setting the cross sec- 
tion to be aei — Tvb^ In A, where the impact parameter is 
b = Ze'^/Ep, Z is the ion's charge, and In A is the Coulomb 
logarithm, wh ich is of order unity in degenerate plasmas 
(|Clavtonlll98"i ). we find 
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(PC); 



(1) 



where the mass density equals lO^pe g cm ^ and fj,e is 
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the electron mean molecular weight. The factor of [1 + 
(pe/Me)^''"^]"^ in equation ((1} accounts for the extension to 
relativistic degeneracy. Equation ([ij shows that degenerate 
plasmas are highly collisional unless the magnetic field ex- 
ceeds > 10^ G (interiors of WDs; pe ~ 1, Me = 2) or > 10^* G 
(interiors of NSs; pa ~ 10*~^, — 10). Hence, magnetically 
driven buoyancy instabilities, if they exist in weakly magne- 
tized NSs or WDs, operate in the collisional limit studied in 
this paper. 

The remainder of this paper is organized as follows. In 
Sj2l we present a linear calculation of buoyancy instabilities 
in weakly magnetized plasmas allowing for arbitrary coUi- 
sionality. We present order of magnitude estimates of the 
growth rates and conditions for instability in the collisional 
limit. We then discuss the application of these instabilities to 
NSs and WDs in 33] We summarize and discuss our results 



2 LINEAR CALCULATION 

We consider a magnetized fluid with an arbitrary equation 
of state and arbitrary coUisionality. The basic equations gov- 
erning the plasma are the continuity equation, 

dp 



dt 



+ V • {pv) = 0, 



(2) 



pg, 



the momentum equation, 

+ p(v ■'V)v = —(V X B) X B - VP - 

at 47r 

the induction equation, 
dB 
dt 

and the energy equation, taken here to be an equation for 
the entropy of the fluid. 



V X (v X B) 



(3) 



(4) 



1 



In Pp 



-ri 



(5) 



where p is the density, v is the fluid velocity, B is the mag- 
netic field, P is the pressure, g is the local gravitational ac- 
celeration, T is the temperature, Fi = {d\nP/ d\n p) ^ and 
Fa — 1 = [d\nT / d\n p) ^ axe adiabatic indices, and s is the 
specific entropy. The heat flux, Q, is due to radiative diffu- 
sion and conduction and is given by 



Q 



-XradVT - XI 



b{b-V))T-xi\h{b-v) T, (6) 



where b = B/B is the unit vector in the direction of B and 
-B = It is helpful to rewrite equation ((Gjl as 



Q 



(Xrad + Xcond) VT - x\\b (b ■ V ) T, 



(7) 



where Xcond ~ X-L is the "isotropic" part of the electron 
conductivity and x'\\ = X|| ~ X-L is then the enhancement 
of the elec t ron co nductivity along the magnetic field. From 
iBraginskiil l| 19651 ). 



X 



X 



and X|| = X, 



(8) 



1 + 

where x is the nonmagnetic thermal conductivity. For 
(^gT 3> 1, the plasma has low coUisionality, x± ~* 0, and 
x'l — X|| = X- This is typically the case for interstellar and 
intergalactic plasmas and stellar coronae. This anisotropic 



conduction gives rise to the buoyancy instabilities described 
in the Introduction. For LjgT <C 1, on the other hand, the 
plasma is highly collisional. This is typical of stellar interi- 
ors, WDs, and weakly magnetized NSs. In this limit, electron 
conduction is largely isotropic with only a small anisotropic 
correction: x± — X and x'\\ — {'^g'^Tx- 

Finally, we assume that there may be composition gra- 
dients in the plasma, but that the timescale for the compo- 
sition to change is long compared to the other timescales of 
interest Q Thus there is no Lagrangian change in the com- 
position: 



-I- V ■ Vpe = 0. 



(9) 



dt 

For concreteness, we have focused on background gradients 
in the electron mean molecular weight /ie, but our calcula- 
tion can easily be generalized to allow for other composition 
gradients (e.g., in the ion abundance). 

We now perform a standard WKBJ analysis on equa- 
tions ([2])- ([9]). We linearize eq uations (|2| - (l9ll using the 
Boussinesq approximation as in iBalbud jloOOl ) and Q08. We 
assume perturbations of the form exp {at — ik ■ x) and that 
gravity is in the vertical direction. Under these assumptions 
the perturbed continuity equation is given by 



fc ■ (5w = 0, 



(10) 



which is the standard incompressibility condition; the per- 
turbed momentum equation is 

.,B-5B .,„B k .SP SpVP 

aSv ^ -ik hiSB ik h — , (11) 

4:Tvp Atvp p P P 

the perturbed induction equation is 

B k 



SB = iSv- 



and the perturbed heat equation is 



(12) 



( -Fia^ + Sv ■ V In Pp-^ A =-ik-5Q (13) 
r3 - 1 V P 1 

where 

= -ixi^okST - <5xcondVr - Sx'iib (b ■ vr) 

-x'l [Sb (b ■ Vr) +b(^Sb- Vt)] - ix'\\b (^k ■ bST^ , (14) 

and Xiso = Xrad + Xcond is the sum of the perpendicular 
electron and radiative conductivities. In the perturbed heat 
equation, we have taken 5P ~ as our implementation of 
the Boussinesq approximation Note that equation (|14p in- 

^ In [O we focus on the application of the collisional MTI to 
the p ^ 10*~® g cm"'' atmospheres of NSs where composition- 
changing reactions are negligible (excluding episodic fusion due 
to accretion). 

^ Taking the dot product of k with equation jlTJl, it is easy to 
show that <5(P -|- /8tt)/P <C Sp/p for short wavelength modes; 
this suggests that we should take 8{P + B^/Stt) in the en- 
ergy equation as our implementation of the Boussinesq approxi- 
mation. However, this introduces spurious growth and damping 
terms into the dispersion relation, as we confirmed by taking the 
incompressible limit, i.e., the high /3 = P/{B^/8tt) limit, of the 
full non-Boussinesq perturbation equations. Thus we implement 
the Boussinesq approximation by taking &P/P <C 5p/p in the 
energy equation. 
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dudes the perturbed conductivities ^Xcond, and We ig- 
nore perturbations in the conductivities proportional to ST 
since these are higher order in 1/kH, where H is the scale 
height. However, this leaves perturbations with respect to B 



-1 2b ■ SB 

B ' 



(15) 



and 5x|| — —(^Xcond- Using equation p5|) in the perturbed 
heat flux, we find 



ikSQ 



Xisofc^ + x'l (b ■ k 



5T + 



i2 



5h\(h- VT ) + ( fc ■ S ) (Sb-VT 



SB ■ b 



1 + Kr)2 B 



[fe- vr-(fe-S)(S-VT) 



(16) 



We calculate ST for any given equation of state, P{p,T, ^l^), 
using the Boussinesq approximation [SP = 0) via 



SP ^0 



dP 



dp / T, 
This implies 
"ainT' 



r f9P 

'^noT 



ST _ 

T ^ 

where 
ainT 



d\np 



Sp 



+ 



aiuT 



91np 
ainT 



d In 



d In pe 

^{dP/dph^ 

T {dP/dT)p,^^ ' 

Pe {dP/dpe)T,p 

' T (9P/9r)p,^, ■ 



ST+ 



Spe 
P,p A*!! 



dP 

dpe 



Spe.{n) 



(18) 

(19) 
(20) 



The ST perturbation has two components: one due to the 
density perturbation and one due to the /le perturbation. To 
calculate Spe, we use the perturbed form of equation (|9]): 



aSpe = Sv ■ V/ie 



(21) 



The temperature perturbation in equation (|18|l induced 
by the density perturbation depends on the thermody- 
namics of the plasma. For instance, equation (I19|) gives 
[dlriT / d\n p]p ^ — —1 for an ideal gas and Spe = 0. On 
the other hand, one can show that for a plasma of non- 
degenerate ions and degenerate electrons, in which the elec- 
trons dominate the pressure (e.g., in a WD or in the outer 
parts of a NS), [dlnT /din p]p = —Pe/riikBT, where Ui is 
the number density of ions. The factor of ZEp/ksT differ- 
ence between the ideal gas and the degenerate plasma arises 
from the fact that the entropy of the degenerate plasma is 
contained in th e ions, which do not contri bute significantly 
to the pressure ijShapiro fc Teukolskvll98^ . Finally, we note 
that at very high densities, when the "ions" (protons and 
neutrons) themselves are degenerate (e.g., in the core of a 
NS), the heat capacity is reduced by an additional factor of 
kBT/Ep^ion, where -Ep,ion is the ion Fermi energy. 

We now combine the perturbed induction equation (eq. 
[12]). the perturbed momentum equation (eq. [H]), and the 
incompressibility condition (eq. |10) ) to find; 



-I- (fc ■ VAf] Sv = -a l^g - ^j-^ 



(22) 



Combining equations JTS}, ([TeJ, ([T8}, ([2TJ, and (dH we ar- 
rive at the final form of the dispersion relation: 



-2 Ar2^ 1 

-aa -oN — + 



dlnT 



dlnp 
dlnT 



-'^cond.ll'?- 



dz 



l^cond, 



+ 

-ki, 

ainT 



91nr 



d In pe 
1 - 2Qbl 
k'z bxk. 



dlnpe fcj 



P.p 



+ 29 



dz k^ 
kxhzkz 



fc2 



1 + (w^r) 



dz 



k-i, k b 



nd,||) 

T~ 

-Prh 
P 

- = (23) 



where k± is the component of the k-vector perpendicular to 
g., i.e. k'i ^ kl + k^, 

1 ainP 91np^ 



TV 



Fi dz 



dz 



and 



is the Brunt- Vaisala frequency, O = 1 — (1 4- [(^gT"]^ 

= + {k ■ va)^ ■ We have also defined the characteristic 
thermal frequencies as 



'^cond.ll — 



Fa - 1 



T 



rxii 



k-bY 



F3-l""PTh 

where Pxh is the thermal component of the total pressure. 

The di spersion relation in equation (|23|l generalizes the 
analyses of lBalbud l|2000l .[ 200 ll ) and Q08 to the case of arbi- 
trary coUisionality and degeneracy. We recover Q08's disper- 
sion relation by taking the coUisionless limit, i.e., cOgT S> 1, 
so that — > 1, and by specializing to an ideal gas equation 
of state and no composition gradients. 

Motivated by the application to NSs and WDs, we now 
assume the following hierarchy of timescales: tdyn ^ (fc • 
va)~^ <S ^^and || ' ^^'^ cousidcr the low frequency (small a) 
expansion of the dispersion relation in equation (|23|l . The 
unstable root of the dispersion relation is then given by 



^ {ijJiso + <JJcond,||) [^st ~ ^^dos] 

^ mk'i/k^ + {k-VAy ' 
where we have defined the stabilizing term as 



(24) 



dlnT 



dlnp 



+ 9 



dlnpe fcj 



dz fc2 



aiuT 



c*ln Pe 



and the (possibly) destabilizing term as 



„2 _ '^cond.ll 
i^dcs = 



UJiBo + ti'cond, 



^ (1 



2ebt) +29 



dlnT PTh 

bz kzbxkx 



fc2 



2(b-fc)~- 



h k 

^^^^ fc2 



k^ 



(26) 



For iV^ > (i.e., a Schwarzschild stable plasma), equation 
(12411 gives the instability criterion: 

(27) 



- fides > 0. 

Note that our definition of Q.st is such that rj^^ 
i) [dlnT/dlnp\p^ 



< since 

< and ii) din pe/dz < for con- 
vective stability (in the sense of the Schwarzschild crite- 
rion, A*'^ > 0). Equation (I27|) is also the general condition 
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^ 

Hurwitz criterion as in iBalbusI { 2000|). Physically, equation 
l|27p is the condition that the destabihzing thermal effects of 
anisotropic conduction in the coUisional limit exceed the sta- 
bilizing effects of magnetic tension and a stable jie gradient. 
The /ie gradient is particularly stabilizing in a degenerate 
plasma. First, it increases the magnitude of A'^^, thus sup- 
pressing the growth rate (eq. [24]). Even more importantly, 
the effect of a stabihzing /^e gradient is at least a factor of 
ZEp/kBT larger than a comparable destabilizing thermal 
gradient. Thus in a degenerate plasma, instability can only 
arise when the /^e gradient is essentially zero. For now, we 
assume that this is the case. We will return to this question 
in 31 

Focusing on the case of zero /i^ gradient, we find that 
the condition for instability is 



(fe • va) - 



aiuT 



c*lnp 



< 0. 



(28) 



For a coUisional plasma, i.e., ujgT <g 1, the instability crite- 
rion becomes, at the order of magnitude level. 



(fe ■ '"a) < ^^gT- 



2Xcond(fc ■ bf dluT 



dz 



dlnT 



dlnp 



Taking {k ■ b)^ ^ k^ , we can rewrite equation (|29|l as 



c; tr^2 ^ 47rpe 2 Xcond „2 



91nr 



dz 



dlnT 



d hi p 



(29) 



(30) 



where H is a typical vertical length-scale in the system, i.e., 
the pressure scale height. Note that the instability criterion 
in equation (|30p is independent of B because the that 
arises from v\ is canceled by the factor of from cOg. The 
order of magnitude growth time for the fastest growing mode 
can now be determined from equation (|24[l : 



■^cond.ll 



tTh [(iJgrKkH)]- 



(31) 



where the largest allowed value of kH is given by equation 
H30p and fxh ~ H'^ / Xcond is the local thermal time. Unlike 
the instability criterion, the growth time does depend on B, 
with tgr oc B~^ for a coUisional plasma. 

Before applying our results to NSs and WDs, we com- 
ment on some of the physics of the pure MTI [bx ~ 1) 
and HBI {bz — 1) instabilities in the coUisional limit. The 
primary difference in our dispersion relation relative to pre- 
vious works is the fact that the conductivity depends on B, 
which introduces terms cx Sx- For the pure HBI case, how- 
ever, where bz — 1 and dT/dz > 0, we find that all terms 
associated with 5x exactly cancel. Thus the basic physics of 
the HBI in Q08 is recovered in the coUisional limit, albeit 
with modifications due to isotropic heat transport associ- 
ated with ojiso. On the other hand, the pure MTI case, in 
which bx = 1, has somewhat different physics in the coUi- 
sional limit. For a low coUisionality plasma, the growth rate 
in the pure MTI limit, for our ordering of timescales, is given 
by 

g\d\nT/d\nz\ 



(32) 



By contrast, the growth rate in a coUisional plasma is 



'^cond.ll ■ 



g\d\nT/d\nz\ 



1 + 



2fc 



(33) 



Equation (|33} shows that the driving of the MTI is enhanced 
in a coUisional plasma relative to what a simple extrapola- 
tion of the low-collisionality result might suggest. Namely, 
we find that including the fact that the magnitude of the 
anisotropy of the heat flux on B in our calculation increases 
the growth rate relative to a calculation that considers only 
a fixed degree of anisotropy (independent of B). Note, how- 
ever, that, despite the appearance of equation (|33p . we re- 
quire fc ■ b 7^ for instability because ii^cond,|| c< {k ■ for 
the pure MTI limit, b^ = 1 and thus growth still requires 
fcx 7^ in a coUisional plasma. Thus the additional driving 
does not change the conditions under which there is growth 
(e.g., dT/dz < 0, weak field, and fcj_ 7^ 0), only the growth 
rate. 

The enhanced driving of the MTI in a coUisional plasma 
arises because of the dependence of the conductivity on B; 
this is analogou s to a k mechanism in stellar oscillation the- 
ory l|Coxlll983 ). To illustrate this point, we study the per- 
turbed heat fluxes for the pure MTI case, which are given 
by 

dT 

-iXisok^ST ~ ix\\kx5T - k^x'w^^^, (34) 



SQz 



-iXisokzST + 2x'h fezC; 



dT 
dz ' 



(35) 



where ^ — 5v/a. In the low-collisionality case studied by 
Q08, the SQ.^ term is the same (except for the inclusion 
of a perturbed isotropic heat flux), while 5Q^ — 0. The 
first term on the RHS of equation (|35p is the perturbed 
isotropic diffusion of heat due to the perturbed tempera- 
ture. The second term is the key new physics: it represents 
the change in the vertical heat flux due to the perturbed 
isotropic conductivity that arises from the perturbed mag- 
netic field, i.e., Sxcoud oc SB. The physical interpretation 
is that as a fiuid element in an MTI unstable situation is 
perturbed, the strength of the magnetic field increases. This 
makes conduction more anisotropic, which increases the con- 
ductive heating/cooling of the fiuid. This in turn leads to a 
buoyant response followed by a greater bending of field lines, 
greater anisotropic conductivity, and yet more heating. This 
additional driving due to the dependence of the conductivity 
on B is the origin of the modest enhancement in the growth 
of the MTI in a coUisional plasma. 



2.1 The Low-collisionality Limit 

For sufficiently strong magneti c fields, a ; qT > 1, an d the low- 
collisi onality limit studied by iBalbusI (|200ll ) and lQuataer"3 
(|2008l ) is appropriate, rather than the coUisional limit high- 
lighted here. For the applications of interest in this paper, 
idyn tTh, in which case the order of magnitude instabil- 
ity criterion in the low-collisionality limit can be determined 
from equations (I25|l - (|27|) with 6 ~* 1. This gives 



{kHf 



(fe ■ bf gH"^ 



dlnT 



dz 



dlnT 



dlnp 



P,Me 



(36) 



Since equation (|36p implies that kmax oc B^^ in the low- 
collisionality limit, decreasing the magnetic fleld strength 
increases the number of unstable modes that can flt in the 
system. This trend continues until the magnetic fleld is suf- 
ficiently weak that ujgT < 1 and the coUisional limit, rather 
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than the coUisionless limit, is appropriate. In the collisional 
limit, the magnetic field dependence of the thermal conduc- 
tivity cancels that of the Alfven speed and kmax is inde- 
pendent of B (eq. [30]). The collisional limit thus allows the 
largest range of wavelengths to be unstable. 

The order of magnitude growth time in the low- 
collisionality limit is given by 



(37) 



The bound on k from equation (|36|l implies that the fastest 
growing mode in the low-coUisionality limit has a growth 
time tgr oc . By contrast, in the collisional limit, the 
fastest growing mode has tgr oc B~^ (eq. [31)). Thus the 
growth time of the MTI and/or HBI is minimized when 
LOgT ^ 1, i.e., at the transition between the collisional and 
low-coUisionality limits. 



3 APPLICATIONS TO COMPACT OBJECTS 

The interiors of main sequence stars have tjgr <§; 1 and 
are thus a possible site for the application of the collisional 
MTI and HBI instabilities. However, it is straightforward 
to show that the order of magnitude instability condition 
(eq. [30]) becomes kH < 10"'^ in the solar interior; thus no 
unstable modes can fit in the star. A more promising site 
for the application of these instabilities is to WDs and NSs, 
since thermal conduction dominates over radiative diffusion 
throughout much of the star. Provided that the instability 
criterion (eq. [3D]) can be satisfied and the growth time (eq. 
|31) ) is sufficiently fast, the collisional MTI and HBI may be 
able to reorient small scale magnetic fields in NSs and WDs 
or provide a mechanism by which the thermal energy can be 
tapped to amplify the local magnetic field. 

3.1 The Cooling Cores of WDs and NSs 

WDs and NSs are born hot, with significant temperature 
gradients left over from the post-main sequence evolution 
of their progenitors. After about a thermal time, the core 
becomes isothermal, with the temperature of the core then 
decreasing in time as energy is lost through the thermally 
insulating surface layers. A natural location to study the 
collisional MTI and HBI is in the cooling core of a WD 
or NS before it becomes isothermal[3 The condition for in- 
stability given by equation (I30|) can be utilized to estimate 
the most unstable wavelength; we assume g « GM/E? and 
dlnT/dz ~ for simplicity. Us i ng th e electron-ion colli- 
sion time from ' Yakovlev fc Urpir] (|l98Cf ) - which is consis- 
tent with equation ([1]) - we find that the unstable modes in 
the core of a nonrelativistic degenerate WD must satisfy 



kR < 13ZK"^/^lnA"^Afn{fi?;7'r«''Ve'^^ 



(38) 



where T = 10 K is the initial temperature of the core of 
the WD, R = 10^i?9 cm is its radius, and M = O.GMo.sMq 
is its mass. 

^ Neutrino cooling at the center of the progenitor star typically 
leads to a temperature maximum at a non-zero radius in newly 
formed WDs, so that the core has both dT/dr < and dT/dr > 0. 
The WD is thus in principle susceptible to both the collisional 
MTI and HBI. The same is true for cooling NSs. 



In the core of a NS, neutron and proton degeneracy 
pressure play an important role. The heat capacity of the 
ions is thus suppressed by a factor of ~ {ksT / E-p Aonf' , which 
reduces the effectiveness of the MTI and HBI relative to 
WDs. We find that modes are unstable in the core of a NS 
provided that 



kR < 0.2Rl^^M-l^''^il%TsZ-''^^ In A" 



(39) 



This shows that the cores of weakly magnetized NSs are 
unlikely to be unstable to the MTI or HBI because the un- 
stable modes have wavelengths larger than the size of the 
star. This estimate did not include the isotropic conduc- 
tivity due to n eutrons in the co r e of the NS, which may 
be substantial (jBaiko et al.ll200ll : iGnedin et al.ll2001^ : this 
would further inhibit the instability. Strongly magnetized 
NSs with B > 10^* G (i.e., magnetars) have ujgT > 1 (eq. 
[T]) and thus the collisional limit in equation (|39|l does not 
apply. However, as shown in H2. li the collisional limit allows 
the largest range of wavelengths to be unstable. Thus the 
MTI and HBI are even less likely to be important in the 
cores of newly formed magnetars. 

Equation (|38p indicates that a cooling WD is unstable 
to the collisional MTI and HBI for long wavelengths > 0.1 R. 
However, the instability in the core of the WD is driven by 
the initial temperature gradient, which vanishes after one 
thermal time, tTh ~ '^^nd' a-fter which the core becomes 
isothermal. Thus growth is only important if tgr < txh- Us- 
ing the estimate of tgr in equation (|3ip , we find the minimum 
magnetic field required for growth: 



B>5x 10^ Z^^' \nK^M^y^R<,T: 



1/2,-1/3 p, 
Me, 2 



(40) 



The strongest magnetic fields measured in WDs are ~ 10^ G 
(Schmidt et al 2003; VanlandiuRiiam ct al. 2005). The large 
magnetic field required for significant growth in the core of 
a WD (eq. [40]) impfies that the MTI and HBI are only 
potentially important for modifying the magnetic fields of 
the most strongly magnetized WDs. 



3.2 The Collisional MTI in the Oceans of 
Neutron Stars 

Even after the core of a WD/NS has become isothermal, a 
temperature gradient can persist in the atmosphere of the 
star. The MTI can operate in this outer envelope. We focus 
on the application to NSs because our estimates indicate 
that the growth times exceed the Hubble time in the atmo- 
spheres of WDs. 

To study the stability properties of NS envelopes we 
first construct hydrostatic constant-flux equilibrium models 
of a NS in plane parallel geometry. We focus on weakly mag- 
netized, highly ionized plasmas and assume that the radia- 
tive opacity is given by a combination of free-free absorp- 
tion and Thomson scattering, while t he conductive opac- 
ity is given by elec tron-ion s cattering (jVentura fc PotekhinI 



20011: iBrown et al.l [2002: Ch ang fc Bildstenll2003l '). foUowing 
Yakovlev fc UrpinlTl98a )l^ We have also tried more modern 
conductivities such as those of jPotekhin et"al] (|199a l: the 



* We have not included magnetic corrections to the opacity in 
the atmosphere model. 
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Igio P[g cm-3] lgj„ p[g cm 3] 

Figure 1. Left: Order of magnitude growtli time (solid line) for the collisional MTI as a function of density in the envelope of a NS with 
Teff = 10® K and B = 10* G. The minimum growth time (maximum growth rate) is fa 6 yrs and occurs where the opacity transitions 
from radiative to conductive (dotted line). In the collisional limit, tgr oc (eq. 1311 ). Right: The maximum wavevector (minimum 

wavelength; eq. |30j ) relative to the pressure scale height H as a function of p for the same model. In the collisional limit, fcmax is 
independent of B. 



conductivities are consistent to about 30% and the growth 
times to a factor of ~ 2, which is sufficient for our purposes. 

In the left panel of Figure [1] we plot the order of magni- 
tude growth time for the fastest growing mode, tgr (eq. [31j 
and [37]), as a function of density in the atmosphere of a 1.4 
M0, 10 km NS with Tos = 10"^ K for B = 10** G (solid line) 
and B = 10^^ G (dashed line). Recall that the growth time 
in the collisional limit scales as tgr oc B~^, while in the col- 
lisionless limit it scales as tgr oc B^. The minimum growth 
time tgr,min is - 50 yrs (~ 0.1 yrs) for 5 = 10* G (10^^ G) 
and occurs at the "sensitivity strip" (jVentura fc PotekhinI 
I2OOII ) where the radiative envelope transitions to the isother- 
mal interior, i.e., where Xcond/xiso ~ 0.5. This is not surpris- 
ing since at smaller radii (higher p), there is no significant 
temperature gradient to drive the instability, while at larger 
radii (lower p), the opacity is largely radiative. The Coulomb 
logarithm F is < 173 where the growth peaks, so that the 
instability occurs largely in the liquid atmosphere, not the 
crust. 

In the right panel of Figure [TJ we plot the most unstable 
wavevector (eq. [30] and [36]) for B = 10* G and 10^^ G, 
respectively. The most unstable wavelength is larger than 
the scale height by at least a factor of 10 for the collisional 
case {B = 10* G) and by a factor of > 300 for the low 
collisionality case {B = 10^^ G). Unlike in the core of the 
star, this does not completely preclude the MTI from act- 
ing because the horizontal wavelength of the mode can be 
much larger than the vertical wavelength in a thin atmo- 
sphere. However, because the tension force is oc k • B, Fig- 
ure [1] does indicate that the instability is most likely to act 
when the background magnetic field is relatively perpendic- 
ular to gravity, so that nearly horizontal perturbations do 
not significantly bend the magnetic field lines. Moreover, 
in H2.ll we showed that the most unstable wavelength for 
the collisional MTI, which is independent of B, is smaller 
than the most unstable wavelength for the MTI in the low- 



colhsionality limit. Thus the B — 10* G result for kma^H in 
Figure [T] is the shortest wavelength for the MTI in any NS 
atmosphere. Since kuia^H < 0.1 even in this case, magnetic 
tension severely limits the extent to which the MTI can op- 
erate in the atmospheres of NSs. In particular, even if the 
MTI is present, tension will significantly modify its nonlinear 
evolution relative to what one would find for kH 3> 1: it will 
be difficult for the field to become relatively radial, as occurs 
in local simulations of t he MTI in the low-collisionality limit 
ijParrish fc Stonell2007l ). 

In spite of the suppression of the MTI by magnetic ten- 
sion, it is nonetheless useful to study a range of NS atmo- 
spheres, parameterized by their effective temperatures Tcft, 
to see how the properties of the instability change in dif- 
ferent NSs. We show these results in Table 1 for B = 10* 
G, i.e., in the collisional limit. We find that the character- 
istic instability time is always short (< 10^ yrs) and that 
the region of maximum growth is always a liquid (F < 173). 
For comparison to our estimated growth times, the typical 
thermal time of a NS at ~ 10*' K is « 10® yrs[3 Thus, at a 
given phase in the evolution of a cooling NS, it is unstable to 
the MTI on a timescale that is very short compared to the 
thermal time on which the properties of the NS evolve. The 
biggest caveat is that, as discussed above, the most unstable 
wavelength at the location where the growth peaks (eq. [30| ) 
is kraaxH ~ 0.01 — 0.1 (Table 1); fcmax -ff is the largest, and 
thus growth is the most hkely, for younger, hotter NSs. 



4 DISCUSSION 

We have studied buoyancy instabilities due to anisotropic 
heat conduction in collisional, degenerate plasmas, i.e., when 

^ This is primarily set by neutrino cooling. Once neutrino cooling 
ceases, the thermal time is > 10^ yrs. 
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Table 1. Table of MTI properties in NS atmo- 
spheres with different effective temperatures 
T(,fl . The columns include the effective temper- 
ature, minimum growth time for the MTI (eq. 
|31|). the density at which the fastest growth 
occurs, the most unstable wavelength in units 
of the local pressure scale height (eq. 1301 ). and 
the local Coulomb parameter. Units are K for 
Tes, yrs for tgr,min and g cm-3 for p^i„. All 
models are for a 1.4 Mq, 10 km NS with B = 
10* G. t„r oc B-^ in the collisional limit. 
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0.089 


0.00653 
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0.0172 


44.0 


6 


96.9 
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3.15 
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10 


47.9 


5.53 
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26.3 
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38.5 
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23.7 


14 


32.8 


12.8 
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21.7 


16 


29.0 


17.3 


0.113 


20.0 


18 


26.6 


22.8 


0.132 


18.6 


20 


24.9 


29.1 


0.151 


17.5 



the electron collision frequency is large compared to the elec- 
tron cyclotron frequency. Although heat conduction is nearly 
isotropic in this limit, i.e., it is nearly independent of B, the 
small residual anisotropy d rives convective instabilities anal- 
ogous to the MTI ((Balbt3l200d ') and HBI (Q08) that have 
been studied previously in low-coUisionality plasmas. The 
physics of the two instabilities is essentially the same in a 
collisional plasma, although there is additional driving of 
the MTI due to the magnetic-field dependence of the ther- 
mal conductivity (see the discussion near eq. [33j 1 . 

In a low-coUisionality plasma, the condition for the 
MTI and/or HBI to grow depends on the magnetic field 
strength B, with magnetic tension suppressing the instabil- 
ity for larger B ( i]2.1[) . The growth time of the instability 
also depends on B, with the fastest growing mode having 
tgr oc B^ . By contrast, in a collisional plasma, the com- 
petition between tension and the destabilizing effects of a 
magnetic-field dependent thermal conductivity leads to an 
instability criterion that is independent of B (eq. [30]); the 
growth time in the collisional limit depends on B, however, 
with tgr oc B~^ . Magnetic tension thus has the smallest ef- 
fect on the MTI and HBI for marginally collisional plasmas 
in which the electron collision frequency is comparable to 
the electron cyclotron frequency. 

The collisional MTI and HBI can in principle operate 
in the cores of young WDs and NSs and in the outer en- 
velopes of NSs (33|. However, magnetic tension and the low 
specific heat of a degenerate plasma severely limit the im- 
portance of the MTI and HBI in these environments. The 
MTI and HBI can undergo several e-foldings in the core of a 
young, high magnetic field WD (> 5 x 10* G) , before thermal 
conduction wipes out the temperature gradient that drives 
the instability in the first place (eq. [40]); however, only 
a small fraction of WDs h ave magnetic fields this strong 
l|Vanlandingham et al.|[2"005l ). We find that there is very un- 
likely to be analogous growth of the MTI or HBI in the cores 



of young NSs (eq. [39]). However, the MTI may be present 
in the liquid oceans of NSs (see Fig. [TJ , where a finite tem- 
perature gradient persists even after the core has become 
isothermal. The most promising candidates are young (hot), 
weakly magnetized (B < 10^ G) NSs. Even in this case, 
however, the fastest growing modes have fcmax-ff ~ 0.1 and 
are thus restricted to very long wavelengths (Table 1). For 
more strongly magnetized NSs, in which conduction is more 
anisotropic, growth in the outer atmosphere is less likely be- 
cause of the increased stabilization by magnetic tension. 

Given these results, one interesting possibility is that 
the collisional MTI could contribute to amplifying and/or 
modifying the magnetic fields in the atmospheres of weakly 
magnetized NSs such as millisecond pulsars or NSs in low- 
mass x-ray binaries. In these systems, it is believed that an 
initially strong birth magnetic field either decayed away or 
was buried by accretion, leaving the NS with a weak field 
that is potentially unstable to the MTI. 

However, it is unclear whether the collisional MTI will 
grow significantly in this context and, if it does, how it 
will saturate. In weak-field simulations of the MTI in low- 
coUisionality plasmas, the magnetic field is amplified by a 
factor of ~ 10 — 30 and the instability saturates by re- 
orienting the magnetic field to be l argely in the direction 
of th e background gravitational field jParrish fc Stonj|2005l . 
|2003). In the oceans of NSs, however, the instability is re- 
stricted to very long wavelengths, larger than the local pres- 
sure scale height of the atmosphere (Fig. [T] and Table 1). 
Thus the instability is most likely to act in regions where 
the local magnetic field is relatively perpendicular to grav- 
ity, so that modes with long horizontal wavelengths do not 
significantly bend the magnetic field. Given these stringent 
requirements for growth, it is difficult to see how the MTI 
can significantly change the structure of the NSs magnetic 
field. Nonlinear simulations are required to assess this in de- 
tail and to determine if the additional driving due to the 
magnetic field-dependent thermal conductivity modifies the 
nonlinear evolution of the MTI in collisional plasmas. 

A background gradient in the composition, in particu- 
lar in the electron mean molecular weight /le, can also act 
to stabilize the MTI and does so particularly effectively in 
a degenerate plasma (see the discussion below eq. [27] )• In 
the part of a NSs atmosphere where the instability is the 
most powerful, however, we do not expect a significant 
gradient. The diffusion time through thi s region is incredi- 
bly f ast, ~ 0.1 - 10 s at p ~ lO'^ gcm"^ (|Chang fc BildstenI 
12004 ). Thus a significant /Xe gradient is only likely if the 
growth happens to occur near a compositional discontinuity 
in the NS. This is not impossible but is unlikely to generi- 
cally be the case. 
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